arXiv: 1501.05737v2 [math.OC] 19 Apr 2016 


STABILIZATION OF PARABOLIC SEMILINEAR EQUATIONS 


lONUT MUNTEANU 


Abstract. We design here a finite-dimensional feedback stabilizing Dirichlet 
boundary controller for the equilibrium solutions to parabolic equations. These 
results extend that ones in which provide a feedback controller expressed 
in terms of the eigenfunctions (pj corresponding to the unstable eigenvalues 
operator corresponding to the linearized equation. In the 
stabilizability result is conditioned by the require of linear independence of 
part of the boundary where control acts. In this work, we 

design a similar control as in and show that it assures the stability of the 
system. This time, we drop the require of linear independence and any other 
additional hypothesis. Some examples are provided in order to illustrate the 
acquired results. More exactly, boundary stabilization of the heat equation 
and the Fitzhugh-Nagumo equation. 
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1. INTRODUCTION 

We consider the following parabolic boundary stabilization problem in a bounded 
domain C with a smooth boundary dfl^ split in two, i.e., dQ = Fi ur 2 , such 
that the Lebesgue measure cr(ri) ^ 0 : 

Pt{t, x) - Ap{t, x) + f{x,p) =0, t > 0, X efi, 
p{t, x) = u{t, x), t > 0, X € Fi, 


^p{t,x) =0, t > 0, X G r2, 


( 1 . 1 ) 


dv 


p{0,x) = Po{x), X Gfl. 


Here, v stands for the unit outward normal on the boundary 512; po G L^(f2) is given 
(the initial data); / is a nonlinear function for which, depending on the context, we 
choose from the following two hypotheses 

(/i) fJpGCinxR); 

(/ 2 ) f,fp,fpp & C'(f2 X R), and there exist Ci > 0, q GN, > 0, i = 
when d = 1, 2, and 0 < < 1, i = 1,..., q, when d = 3, such that 

\fppix,p)\<Ci (^|p|“^ + l) , VpeR. 


Here, fp, fpp denote the partial derivatives of / with respect to its second argument. 
On the part of the boundary Fi it is applied the control u, while F 2 is insulated. 
Let pe G (7^(12) be an equilibrium solution to (|1.1|) . that is, pe satisfies 

d 

-Ape + f{x,Pe) = 0 in 12; pe = 0 on Fi, —pe = 0 on F 2 . 

av 
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Then, defining the fluctuation variable w := p — pe, equation (HID can be rewritten 
as 


Wt{t,x) - Aw{t,x) + f{x,W +Pe) - f{x,Pe) =0, t > 0, X € O, 
w{t, x) =v := u{t, x) — Pe{x), t > 0, a: G Ti, 

■^w{t,x) = 0, t > 0, cc G r2, 

r(;(0, x) = Woix) := Po{x) - Pe{x), x € VI. 


Our main concern here is the design of a boundary feedback controller u that 
stabilizes exponentially the equilibrium state Pe in (HD, or, equivalently, the zero 
solution in m- 

The problem of boundary feedback stabilization of this kind of parabolic-like 
equations was first solved in the pioneering work [14) . Then, several others meth¬ 
ods were proposed for deriving new types of controls, like backstepping approach, 
see maiiiig, or the proportional type controllers designed in [I m 113] . The last 
ones provide a simple form feedback stabilizer expressed in terms of the unstable 
eigenfunctions of the operator given by the linearized equation around the equilib¬ 
rium Pe (that is the operator A below), which is easily implementable in practice. 
However, the stabilizing procedure is conditioned by the require that the system 
linearly independent in where, (j)i, i = l,...,iV, are the cor¬ 

responding first N eigenfunctions. In the present work, we extend these results in 
the sense that we drop the hypothesis of linear independence (which, in fact, seems 
to be quite restrictive, as shown in m) . More precisely, we design a similar type of 
control as that one in [1] , and show that it assures the stability of the steady-state 
Pe, without needing the linear independence assumption. Roughly saying, the idea 

f 'V 

is as follows: in the case when the system \ -§p4’is is not linearly independent, 
the corresponding Gram matrix is singular, and, consequently, the approach in [T] 
fails because the controller is not well-defined anymore. However, we overcome 
this in the next manner: for the beginning, in order to clarify the procedure, we 
work under the assumption that the first fV G N unstable eigenvalues are simple. 
We involve N matrices obtained from the Gramian, multiplied by some diagonal 
matrices, and show that their sum is invertible, even if each one is not. Therefore, 
we may well-define u as the sum of N different controllers obtained from those 
matrices, and arrive to similar stabilization results as in [I]. The case of general 
eigenvalues involves a similar boundary feedback, as the one described above. The 
idea is to artificially make the eigenvalues simple again by slightly changing the 
linear operator. Finally, via a fixed point argument, locally stabilization results 
may be deduced for the fully nonlinear system (HD. 

On this subject, it should be mentioned also the work [6], where controllers, of 
the above type, are designed. Instead of the eigenfunctions system, they show that 
there exists a (possible different) family of functions for which, the corresponding 
proportional controller assures the stability. That result is based on some unique 
continuation property for the involved operators. Also, the work m , on the Navier- 
Stokes equations, provide similar proportional feedbacks, but with the coefficients 
given by Riccati algebraic equations. However, in comparison to those results, here, 
the form of the feedback is given explicitly, and no additional assumption to the 
classical context is needed. 
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2. The design of the feedback law and the main results 

In this section, we detail the design of the new feedback, show how it differs from 
that one in [T], and state the main stabilization results. To this end, let us define 
the operator A : V{A) —>■ as 

Ay := -Ay + a{x)y, Vy S V(A) = |y € : y|ri = 0, ^y\r 2 = o| , 

where a{x) := fy{x,yeix)), x € ft. Notice that A is self-adjoint, besides this, it 
can be shown that A has compact resolvent, therefore, it has a countable set of 
eigenvalues, denoted by (repeated accordingly to their multiplicity) such 

that, given p > 0, there exists only a finite number of eigenvalues with 

Xi < p, i = (we call them the unstable eigenvalues). So, Ai > p, for all 

i > iV-l- 1 (for more details, see [I]). As we see below, the larger p is, the faster the 
exponential decay of the solution is. (But, on the other hand, the larger p is, the 
larger dimension, A^, of the feedback is.) Let us denote by {4>i}'^i the corresponding 
eigenfunctions, that can be chosen in such a way to form an orthonormal basis of 
the space T’(A). 

2.1. The case of simple eigenvalues. For the beginning, in order to make clearer 
our approach, we add to the above context the following hypothesis: 

(iJ) The eigenvalues Ai, * = 1,..., N, are mutually distinct. 

In other words, we assume that the first N eigenvalues are simple. So, we can 
(eventually) rearrange them such that 

Ai < A2 < ... < Ajv. 

Then, in the subsection below, we show how we may drop this assumption, in order 
to obtain stabilization results in the framework of general eigenvalues. 

One of the tricks used here is to ’’lift” the boundary conditions into the equations. 
For this, we introduce the so-called Dirichlet operator as: given a G T^(ri) and 
7 > 0, we denote by D^a := y, the solution to the equation 

N 

-Ay{x) + a{x)y{x) - 2^ Afc {y,(j)k) 4>k{x) + ^y{x) = 0, for x 

k=i (2.1) 

y = a, on Fi, 

|;y = 0, onr2. 

For 7 > 0 large enough equation (12.11) has a solution, defining so the map D.^ £ 
L(L^{ri),H^{ft)). Besides this, we have 

||A’7a||^i < —||Q;||L2(rj), Va G L^(ri). 

(see, e.g., m p-6> line 16]) We choose p < 71 < 72 < ... < 7Ar, constants such 
that equation m is well-posed for each of them, and denote by D.y., i = 1,..., A^, 
the corresponding Dirichlet operator. 

Next, let us denote by B the Gram matrix of the system | in the Hilbert 

space L^(Fi), with the standard scalar product (y,/i)g := /p^ f{x)g{x)da [a being 
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the corresponding Lebesgue measure on the boundary Fi). That is 


B := 




Further, we introduce the matrices 


/ 1 

7k-M 


^Ik 


Ik—M 

0 


0 

0 

1 

lk—>~N 




, k = 


( 7i-Ai 
1 


T ■= 


V 


72 —Ai dv ' 


^|7'^2(a:) 




and 


where 


A = (Bi + B 2 + ... + Bjv) 

Bfe := fc=l,...,Af. 


( 2 . 2 ) 


(2.3) 


(2.4) 


(2.5) 

( 2 . 6 ) 


At this stage, we are able to give the form of the stabilizing feedback, but, before 
this let us take a pause and make some comments on the invertibility of the sum 
_Bi + ... + Bjv, which is, in fact, the key result of this paper. More exactly, since the 
system is not necessarily linearly independent in L^(ri), the gramian B 

may be singular, and, therefore, each Bk,k = may be singular. However, 

by Lemma 15.21 in the Appendix, the sum Bi + ... + B^ is an invertible matrix, 
provided that hypothesis {H) holds true. Therefore, the matrix A is well-defined. 

Now, let us introduce the feedbacks 


Uk{t,x) 


A 


( {p{t) - Pe,(j)l) \ 
{p{t) -pe,(/>2) 


/ \ 


+ Pe, t > 0, a: G Fi, 


\{p{t)-p,,Mj \ J ^ 

(2.7) 

for k = Here (•,•) denotes the standard scalar product in L^(H); while 

(•, •)jY denotes the standard scalar product in R^. 

In the case where the Gram matrix B is nonsingular, one may consider only 
one feedback from the above list, for example ui (since, in this case, B^^ exists), 
take u = ui and show that it assures the stability stated in the main results of the 
paper, Theorems 12.1112.21 below. In this case, in fact, we stumble on the feedback 
designed in [I]. 

When B is not invertible anymore, we take u to be the sum 


U = Ml -I- U2 -I- ... + UN, 
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which, in a condensed form, can be written as 


,= iT A 


\ir \ / IT T ^ ! 

{pit) -Pe, 4 > 2 ) 

\ {pit) -Pe,<l)N) ) 


\ 

1 

1 

y 

1 1 y 


-Pe- 


( 2 . 8 ) 


N 


We claim that the above feedback assures the stability of the steady-state pe in 
(HH). In what follows, we focus to prove this is indeed so. The main ingredient to¬ 
ward the proposed goal is the stabilization of the linearized equation corresponding 
to m, given by. 


yt{t, x) - Ay{t, x) -|- a{x)y{t, x) = 0, t > 0, a; G 11, 
y{t,x) = v(l,x), t > 0, X G Ti, 

■§^y{t,x) =0, t > 0, X G Ta, 

?/(0,x) = yo, X G fl. 


The following results amounts to saying that the feedback u achieves global expo¬ 
nential stability in the linear system (1^ . and local exponential stability in dm. 
More precisely 


Theorem 2.1. Assume that hypothesis (H) holds and f obeys (/i). The feedback u, 
given by h2.t^) . exponentially stabilizes the linearized equation More exactly, 

the solution y to the equation 


yt{t,x) - Ay{t,x) + fy{x,yeix))y{t,x) = 0, t > 0, x G fl, 

/ {y{t),(l)i) \ / 1 \ 

{y{t),(j)2) 


y{t,x) = {T A 


V {yit),M J 

-^y{t,x) =0, t > 0, X G Tz, 


1 

V1 y 


t > 0, X G Ti, 


( 2 . 10 ) 


N 


. y(0,x) =yo, X G fl, 
satisfies the exponential decay 

||y(t)f <Ce-^‘||yof, t > 0, 


( 2 . 11 ) 


for a prescribed p > 0, and a constant C > 0. Here T is introduced in relation 
while A is introduced in relation /fO) . II • II denotes the norm in L'^{fl). 


Theorem 2.2. Assume that hypotheses (i/),(/ 2 ) hold true. Then, the solution to 
the closed-loop nonlinear equation 


pt{t,x) - Ap{t,x) -i- f{x,p) = 0, t > 0, X G 11 
{p{t) -Pe.fil) 

(p{t) -Pe,02) 

V {P{t) - Pe, (I)n) J 
-§^p{t,x) = 0, 1 > 0, X G r2. 


p{t,x) = (t A 


\ 

1 

y 

1 1 y 


-{-Pe{x), t > 0, X G Ti, 


N 


^ p(0,x) =Po, X G H, 

satisfies the exponential decay 

\\p{t) -Pef < Ce~f^*\\po - PeW^, t > 0, 


( 2 . 12 ) 
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for a prescribed /i > 0, and a constant C > 0, provided that \\po — Pe|| is small 
enough. Here T is introduced in relation fS.4D , while A is introduced in relation 

The proofs of these two theorems will be given latter, in Section |4l 


2.2. The case of general eigenvalues. As announced above, we see that, slightly 
perturbing the linear operator A, we are able to show that the present approach 
works for the case of general unstable eigenvalues (i.e., not necessarily simple). 
Indeed, let us assume, for example, that the first unstable eigenvalue has its mul¬ 
tiplicity equal to 2, i.e., Ai = A 2 , and the others unstable eigenvalues are simple 
(other cases can be treated in a similar manner as below). 

This time, the Dirichlet operator is introduced as follows: for any a G L‘^{Ti), 
we define D^a '■= y, where y is solution to 


N 


-Ay{x) + a{x)y{x) - 2 ^ Afc {y, 4>k) (j)k{x) - 6 {y, (fi) (^1 -b jy{x) = 0, x Gfl, 


y = Of on r 1 , 
|^ = 0on r^. 




(2.13) 


for some (5 > 0. We choose p < 71 < 72 := 7 i + < 73 •= 7i + Iv^"- < 7^ 

7 i -b 1, with 7 i large enough, such that, for (5 = ^ equation (12.131) is well-posed for 
each 7 i, i = 1,N, and such that Ai -b <5, A 2 ,..., Xn are mutually distinct. 

Now, our feedback has the form 


u= (T A 


{p{t) -Pe,<t)2) 

V {pit) -Pe,(l}N) J 


\ 

( ' \ 

? 

1 

/ 

11 y 


-Pe- 


N 


T is defined as 


d 


T := 


(fiix) 


^ 71 — ( 5 — Ai dv ^ 

72 -5-Ai (^) 

—^— —Ai{x) 




and 

where 

for 


7 ^ —5—Ai dv 7iV-A2 

A = {Bi + B 2 + ... -b Bn)~^, 
Bk k = l,...,N 


7 Ar-Ajv '§ii^N{x) J 


B := 


and 


\ {^4>n,-§^4>s)q 

/ _ 1 _ 

7fc—(5 —Ai 




\ 


0 

1 

7fc—A 2 

0 


(J;^2,^9i>Ar)o 
{-§d^N , -§p(I)n) Q ) 

\ 


, k = l,...,N. 


7fc —An 


/ 


(2.14) 


(2.15) 

(2.16) 

(2.17) 

(2.18) 


(2.19) 
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By the above definitions, taking into account that Ai + 6, X 2 , ■■■, Xn are mutually 
distinct and 71 , 72 ,are also mutually distinct, one may show that the matrix 
A is well-defined, that is, the sum Bi + B 2 + ■■■ + Bn is indeed invertible (this 
follows similarly as in Lemma l5.II and Lemma 15.21 in the Appendix). 

The main results, in the context of general eigenvalues, are stated below. 


Theorem 2.3. Assume that f satisfies (/i). Then, the feedback u, given by 
exponentially stabilizes the linearized system JO) . More precisely, the solution y 
to the system 


ytit,x) - Ay{t,x) + fy{x,yeix))y{t,x) =0, t > 0, x G fl, 


y{t,x) = (T A 


\ {yit),M J 

-§py{t,x) = 0 , t > 0 , xer 2 , 


1 

V 1 ) 


t > 0, X € Li, 


( 2 . 20 ) 


N 


. 2/(0,a:) = 1 / 0 , X e n, 


satisfies the exponential decay 


\\y{t)r<Ce-^*\\yor, t>0, 


( 2 . 21 ) 


for a prescribed yt > 0, and a constant C > 0. 
i2.15\} . while A is introduced in relation \2.16\) . 


Here T is introduced in relation 
■ II denotes the norm in L^(fX). 


And, concerning the nonlinear equation. 


Theorem 2.4. Let 1 < d < 3 and f obeying assumption (/ 2 ). The solution to the 
closed-loop nonlinear system 

pt{t,x) - Ap{t,x) f{x,p) = 0, t>0, X e LI, 


p(t,x) = It A 


( {p{t) - pe, fi) \ 
{p{t) -Pe,h) 

V {p{t) -Pe,(l>N) J 


—p{t,x) = 0 , t > 0 , X e r 2 . 


( ' \ 
1 

V 1 / 


■TPe{x), t > 0, X e Li, 


N 


dv 


( 2 . 22 ) 


. p(0,x) =Po, X e Ll, 

satisfies the exponential decay 

\\p{t) -Pef < Ce~'"*\\po - Pef , t > 0, 

for a prescribed yi > 0, and a constant C > 0, provided that \\po — Pe|| is small 
enough. Here T is introduced in relation Ii2.15\) . while A is introduced in relation 

IKM- 

The proofs will be given latter, in Section |4l 


3. Examples 

In order to illustrate the acquired results, let us consider some examples arising 
from thermodynamics and biology. These examples are treated briefly, since we do 
not want to oversize the present work. However, more details and also numerical 
simulations can be found in the work [12] on the stabilization of the Fischer’s 
equation by similar proportional feedbacks as above. 
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Example 1. We consider first a problem discussed in [4], and also in [I]. More 
precisely, the stabilization of the heat equation on the rod (0,1) 


Vt - Vxx - Ay, X € (0,1), t > 0, 
yx{t, 0) = 0, y{t, 1) = u{t), t > 0, 


(3.1) 


with the Dirichlet actuation u in a; = 1 and with A a positive constant parameter. 
First of all, we notice that the stabilizing results using the backsteppting method 
in [4] work for arbitrarly level of instability, while in [1] they are applicable un¬ 
der the condition that the parameter A do not exceed the bound (^) . In the 
present paper, Theorem 12.11 above holds true without any a priori condition on 
A. This is indeed so, because, in our case, the operator Ay = —y" — Ay, Vy £ 

V^A) = {y £ 1) : y'(0) = y(l) = 0} has the eigenvalues Xj = ^ -^ 

with the corresponding eigenfunctions (j)j = cos j = 1,2,.... It is clear 

that hypothesis (H) holds true for any N G N. Let us pick some N such that 
Xj > p > 0, yj = N + 1, N + 2,.... Since system dSlD is linear and Theorem 12.II is 
applicable we get that the feedback 


/ 


= ( TA 


/j y{^) cos ^xdx 
/p y{x) cos ^xdx 


\ 


V Jo ^^^ xdx J 


( ' \ 
1 

V 1 


N 


where 


/ - 


T = 


27r 

677 

ll —TT^ 

/l — 

2-77 

677 

^2— 

/2-(37r)2 

277 

677 




/_,xAr_2(2A-]>_ 

v /i-[(2Af-l)7r]^ 


(- 1 ) 


N 2(2Ar-l)7r 
Z2-[(2Af-l)7r]2 

2(2jV-l)7r 


t 2^2iV-iJ7r 


and 
A = 


= 


N 

E 

k=l 




_-3_ 

lh-7r^llh-(37r)^l 


(_l)J^+i(2Af-i) 

\ lh-7r^llh-(2JV-l)^7r^] 


-3 

_9 

14-(37I 


TTTT 


(_1)W+13(2JV-1) 

■lZ.-(37r)^|l4-(2A-l)^^^| 


(_i)W+i(2Ar-i) 

[Zfc-7r^][Zfc-(2iV-l)^7rq 

(-l)^3(2Af-l) 

[/fc-(37r)2][is,-(2Ar-l)27r2] 


(2Ar-l)^ 

14_(2JV-l)^77q^ 


(here we have denoted by Ik ■= 4(7fc -b A), k = 1,..., N), assures the global expo¬ 
nential stability of the null solution of (j3.ip . Here, p < 71 < 72 < •.• < 7 Ar are any 
positive constants. 

Example 2. We consider now the classical Fitzhugh-Nagumo equation that de¬ 
scribes the dynamics of electrical potential across cell membrane (for details, see 

i) 

yt - yxx + y{y - l)(y - a) = 0, 0 < a; < Z, t > 0, . . 

yoo{t,l)=0, y(t, 0 )=u(t), t > 0 , ^ 

where 0 < a < ^. It is known that the equilibrium y^ = a is unstable. The 
linearized operator A is given by 

Aly = -y" - a(I - a)y, Vy £ V{A) = {y £ iF^(0, 1) : y(0) = y'(Z) = O} . 
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f 2 11 ^ TT ^ 

It has the simple eigenvalues Xj = ^ -a(l — a), j = 1,2,..., with the 

corresponding eigenfunctions (j)j = sin x, j = 1,2,.... Since hypothesis (H) 

is full filed for any G N, we chose one such that Xj > p > 0, j = N -\-l, N -\-2,.... 
It is easy to see that f{y) = y{y — l){y — a) obeys the hypothesis of Theorem 12.21 
and so we get that the feedback 


u = (TA 


X 


2r 

fg(y(x) — a) sin ^xdx 

V loiyX) - «) sin '^^^^r^xdx j 

( 


where 



( 1 \ 

? 

1 

/ 

11 / 


N 


A = 


T = 


( 


N 


(2l)V£ 


k=l 


21tt 

QItt 

2(2N-l)l \ 


Zl —TT^ 
21‘K 

/i-(37r)2 

S/tt 

Ii-(2N-1)^'I7^] 

2{2N-l)l 



G —(37r)^ 

/2-(2Af-l)^7i-^] 


21tt 

G/tt 

2(2W-l)i 


1N —TT^ 

/n —(37r)2 

iw-(2Ar-l)27r2] / 


- a)), k 

= l,2,...,iV; 

and 


( 1 


3 

2N-1 

V h-'!r 


l4-7r^J|ifc-(3-7r)^| 

|/fc--7r^J|ii,-(2A7-l)^7r^| 

3 



3(2Ar-l) 


-(371-)^] 

Lfc-(37r)-^ J 

lh-(37r)^Jlifc-(2A''-l)^7r^J 


_ 2Af-l _ _ 3(2Ar-l) _ 

\ lh-7r^llh -{ 2 Af-l)^ir^| [4 - (Stt)^] [ 4 -(2A'-l)^7r^] 

locally stabilizes the solution y^, = a in (lO) . 

Example 3. Finally, we consider the periodic heat equation in ( 0 , 7 r)^ 

yt - Ay - fiy = 0, x G (0, Tr)^, t > 0, 
y(t,xi,0) = u{t), yx.,{t,xi,TT) = 0, Xi G ( 0 , 7 r), 

2 / 3,1 (t, 0 ,X 2 ) = 2 / 3,1 (t, 7 r,X 2 ) =0, X 2 G ( 0 , 7 r). 

In this case the operator 

Ay = -Ay - py, 

for all 


( h 


2N-1 

-(2Ar-l)^7r^ 


(3.3) 


y G V{A) = {?/ G X{{0,ttX : y(xi,0) = 2 / 3 ; 2 (a^i, tt) = 0, 2 /xi( 0 ,X 2 ) = 2 /xi( 7 r,X 2 ) = O} , 
has the eigenvalues 

= -/i, Vfc= (fciife) gN^, 

with the corresponding eigenfunctions 

(j)k = cosfcixi sin ^ X 2 , Vfc G N^. 

Ordering the eigenvalues set as an increasing sequence and redefine them, we have 
Ai = 1.25 — p, = cosxisinix 2 ; A 2 = 3.25 —/i, 4>2 = cosxisin|x 2 ; A 3 = 
4.25 — p, (j )3 = cos 2xi sin ix 2 ; A 4 = 6.25 — p, (/>4 = cos 2xi sin |x 2 ; A 5 = 7.25 — 

/r, (/is = cosxi sin |x 2 ; Ae = 9.25 —/r, (/ig = cos 3xi sin ix 2 ; Xy = 10.25 — p,(j)7 = 
cos 2xi sin |x 2 ; Ag = 11.25—p, (pg = cos 3xi sin |x 2 ; Ag = 13.25—/r, pg = cos Xi sin |x 2 ; 
Aio = 15.25 —p, pio = cos3xi sin |x 2 ; An = 16.25 —p, cpn = cos4xi sin ix 2 ; A 12 = 
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16.25 — /J,, 012 = cos 2xi sin ^X 2 - It is clear that since An = A 12 , hypothesis (H) fails 
to hold. So, we apply this time the result in Theorem l2.3l Thus, the corresponding 
control u of the form (|2.14p assures the stability of the null solution in the system 
(1331) . We will not write it down here explicitly since it involves three 11th order 
square matrices. 


4. Proofs of the results 


Proof of Theorem \2.1[ The proof is based on similar ideas with those ones 
in the proof of [1] Theorem 4.1]. We shall show only the fact that the first N modes 
of the solution y are exponentially decaying. The rest will be omitted because it 
can be deduced by almost identical arguments as in the proof of [U Theorem 4.1]. 
The main difference between the present proof and that one in [I] is that while in 
[I], the stability of the unstable modes is showed for each one independently, here 
we consider the system formed by them, and show that it is stable. That is why, 
the present feedback involves matrices, rather than vectors as it does in [ 1 ]. 

Firstly, recall the Dirichlet operators i = defined after relation 

dllD. In the following, we need to compute the scalar product for all 

b j G {!) •■•jTo this end, we scalarly multiply equation (|2.1I1 . corresponding to 
7 j, by 4>i, to get, via Green’s formula, that 


0 = - / l^y{x)4>i{x)dx-\- / a{x)y{x)(j)i{x)dx + -2\i) / y{x)4>i[x)dx 

J Vi J Vl J Vl 

= [ a{x)-^(l)i{x)da + [ j/(a;)(-A0i(a;) + a(a;)0i(x))cia: + ( 7 , - 2Ai) / y{x)(t>i{x)dx. 
Jfi dv Jq Jq 

(4.1) 


It yields that 


a, (pi'^ — 


1 


7j - Ai 

Next, we introduce the feedbacks 


a, ^0*/ ’ ^ I,---,^- 


(4.2) 


Vk{t,x) = ( ^ 

{y{t),h) 

( 

\ 

V {y{t),(pN) ) 

V 


1 0 

7 fc-Ai dv 

1 d 

Jk->‘2 du 


02 (a;) 


t > 0, a; G Ti, (4.3) 


Ik-^N dv 


N 


for fc = 1, 2 ,..., N] and v = + U 2 + ■•■ + vn- It is easy to see that we have 


= (T A 


( ( 2 /(i),</>i) \ 

\ j 


I ' \ 

1 

V 1 / 


N 


that is exactly the boundary feedback plugged in (12.101) . Here, T is defined in (12.4L 
while A is defined in (1351) . 

For latter purpose, we show that 


/ {D^^Vk.pi) \ 
{D^^Vk.h) 

= —BkA 

( {yif),h) \ 
{y{t),h) 

\ {D^^Vk,4'N) y 


V {y{t),(pN) j 


(4.4) 
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where Bk are introduced in (12.611 above, for k = This is indeed so. We 

have 


'Pi) ( ^ 


/ {y{t),<Pi) \ 
{yit),(p2) 

V {y{t),(pN) j 


/ 7fc-Ai \ 

7fc-A2 

^ 7fc-AN {^Ik-g^pNTpi) ) 


i = 1, ...,N. 


N 


It then follows by relation (14.21) . that 

{^~fk''’kT 'pi) ~ \ ^ 


( {y{P)Ai) \ 

{yit),p2) 


{lk—>'i)('yk—>'i) duP^')0 

hk-->'2)(7k--ki) ( ^2) dvPP) 0 


\ {y{t),pN) J 


i = 1 ,..., 


\ {7k-'^N){7k — \i) ^dvP^ ^ dvPP) 0 ) 


N 


from where we immediately obtain (14.41) . 

Now, let us return to the linear system 

yt{t, x) - Ay{t, x) + a{x)y{t, x) = 0, t > 0, x € ft, 
y{t,x) = v{t,x) = vi{t,x) + ... + VN{t,x), t > 0, cceTi, 
^y{t,x) =0, t > 0, x G Ta, 

1 /(0, x) = yo[x), X €ft. 


(4.5) 


(4.6) 


Let us denote by z{t,x) := y{t,x) — D.y^vi{t,x) — ... — D.yi^V]s[{t,x), t >0, x € ft. 
We claim that the feedbacks Vk, k = 1,..., N, can be expressed only in terms of z, 
as 

1 a 


Vk{t,x) = \i^ 


( {z{t),Pl) \ 

{z{t),p2) 

V {z{t),pN) J 


( 'yk-XidvP^ \ 
lk-^2 

^ -Vi.—Aw dvP^ > 


(4.7) 


ik—y^N dv 

To see this we do the following straightforward computations 


AT 


I A 


( {z{t),Pl) \ 
{z{t),p2) 

V {z{t),(l)N) J 


I 7 fc-Ai duP^ \ 

1 a A 

71c- a2 a,Ap'^ 

V 'vc.-Aw ^pN J 


= 2 ^ 


-\tu 


7k — >'N du 


N 
1 d 


( {y{t),Pi) \ 

{yit),p2) 

V {y{t),pN) J 


1 

a 

7(c-Ai 

av 

1 

a 

lk-^2 

av 

1 

a 


7fc-Aiv av 


pN 


i=l 


(Dj,Vt(t),p2) 

V {DjiVi{t),pN) J 


( 7k->.i a^P^ \ 

V -c. — Aw a^P^ ^ 


N 


(taking into account relation (14.41) 1 


— — M7 + A{Bi + ... + Bn)] A 


( {y{t),Pi) \ 
{y{t),p2) 

V {yit),pN) J 


( 7k->^iauP^ \ 

^ ^ p2 


7k-ik2 du 

7k->'N / 


\ -A 

^ ‘yu — 


N 


= Vk, 


-e- -Q- 
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since A = {Bi + ... + Bn) ^. Moreover, likewise in ()4.4I1 . we have now 

= (4.8) 


(4.9) 


1 \ 



{D^Ak,h) 

= -^Bi,a 


\ {D-iAk,4>N) j 


V {z{t),(j)N) J 


Finally, equation (14.61) may be rewritten in terms of z as follows 

f zt{t, x) + Az{t, x) = R{{z, (l)i) ,..., (z, (I)n)), t>0, X efi, 

I z(0,x) = Zo(x), X e fl, 

where 

C N \ N N 

^ Aj (jij + ^ %D^.Vi. 

i=l / t i,j=i i=l 

(4.10) 

By (14.81) . and using the fact that A is the inverse of the sum oi Bk, k = 1,...,N, 
we see immediately that 

/ \ 


{R,h) 

V {r,(I}n) J 


.. N N N 

= - ^ BkAZt + A ^ AZ - - ^ ikBkAZ 

k^l 


1 11 ^ 

= + AZ — 271-2 + - ^(71 “ lk)BkAZ, 

k=2 


(4.11) 


where we have denoted by Z(t) := 

( {z{t),(l>i) \ 
{zit),(l)2) 

f 

, t > 0; and by A := 


V {z{t),(j)N) ) 

V 


0 A 2 ... 


\ 0 0 ... 

We have now, by (gH), scalarly multiplied by (j)ij * = together with 

(I4.11|) . that 

1 1 1 ^ 

Zt + AZ= -Zt +AZ- - 7 iZ + - ^(71 - ^k)BkAZ, t > 0; Z(0) = Z^. (4.12) 

k^2 


Or, equivalently, 

Zt = - 7 iZ + ^(71 - -ik)BkAZ, t > 0; Z(0) = Z„. 


N 


(4.13) 


k=2 


Recall that Bj, j = 1,...,N, are positive semidefinite symmetric matrices (by the 
definition of Bj,K.y. and the fact that B is a Gram matrix), therefore, {Bjq,q)j^ > 
0, Vg S K^, j = 1,..., N. Consequently, A = [Bi + ... + Bn)~^ is a positive definite 
symmetric matrix. Thus one can define another positive definite symmetric matrix, 
denoted by A 2 , such that A 2 A 2 = A (the square root of A; for details see [7]). Let 
us scalarly multiply equation ()4.13ll by AZ, to get 


l_d 
2 dt 


N 


||A^Z(t)||^ = - 7 i||Aiz(t)||^ + ^(71 - 7 fe) {BkAZit), AZ{t))^ , (4.14) 


0 

0 

Xn 


k=2 
















STABILIZATION OF PARABOLIC SEMILINEAR EQUATIONS 


13 


that leads to 

]^±\\A^Z{t)\\l<-^4A^Z(t)\\l, t>0, 

since 71 — 7 fc< 0 , k = 2, Here || • Hat stands for the euclidean norm in . 

The above relation implies the exponential decay of Z in the WA^ ■ ||Ar-norm, i.e., 

where using the fact that Ai is a positive definite symmetric matrix, we finally 
arrive to 

\\Zit)\\j,<Ce-^^^*\\Zo\\%, t>0, (4.15) 

for some positive constant C. 

The rest of the proof follows by identical arguments as in the proof of [U Theorem 
4.1], therefore it is omitted. However, we notice that, as shown in [1], the remaining 
infinitely many modes will be exponentially stabilized by the feedback v with the 
gain e~P*. Hence, by (14.151) . since 71 and p may be taken arbitrarily large, the 
decaying constant p in (| 2 . 1 ip may be taken arbitrarily large as well. □ 

In the case of general eigenvalues, we have 

Proof of Theorem \2.,'^ The proof is almost identical with that one of Theorem 
o that is why we shall skip most of computational details. 

Recall the Dirichlet operators Dj., i = 1, ...,iV, introduced after relation (12.131) . 
Likewise in (l4T])-(IO) we get 


(Dj.a,(j)i) — ^ and {D^.a,(j)i) — (^ 0 ;, at* )p > 


for all i = 2 ,..., N, j = 1,..., N. 
Now, introduce the feedbacks 


Vk{t,x) = { A 


(4.16) 



( 

(2/(t),</'2) 

5 

^ {y{t)AN) y 

V 


1 d 




'yk->'2 dv 

^ ^: 4 'n{x) j 


, t > 0, a; G Ti, (4.17) 


7fc—•''iV \- / / jy 

for k = 1, and v = ui + U 2 + ... + vn- It is clear that 


It A 

( (y(t),0i) \ 
{y{t),h) 

5 

1 

\ 

V {y(t),4>N) ) 


1 1 / 


N 


which is the boundary feedback plugged in ( 12 . 201 ) . 

Set z{t, x) := y{t, x) — x) — Dy^V 2 {t, x) — ... — Dy^VNit, x). Similarly as 

in (IQ) we have 


Vk{t,x) 




( 


•> 


V 


■yk-S-Xi \ 

lk-A 2 OU 


jk->^N du ' 


N 


k = I,...,iV. 


(4.18) 
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Besides this, as in (|4.8|) . we also have that 


^ {D-fkVk,4i) '' 


( {zit),4i) \ 

{D^kVk,42) 

= -^BkA 

{z{t),42) 

\ {D^Ak,4N) y 


V {z{t),4N) j 


k = l,...,N. 


(4.19) 


In terms of the new variable z, the linearization (j2.20|) may be rewritten as 

zt{t,x) + Az{t,x) = R{{z,4) ,{z,(j) 3 ),..., {z,(j)N)), t>o, X en, 
z{0,x) = Zoix) := Tjo - D^^vi{Q,x) - 1)2(0,a;) - ... - D^^vn{0,x), x G O, 

(4.20) 

where 


N 


N 


R{{z, 4i ). {z, h), ■■■, {z, (I)n)) := - ( X! ^ Aj {D^.V^, 4^) 4j 

N N 

Then, simple computations as in f|4.1ip give 

/ (^.<^ 1 ) \ 

{R, 42 ) ^ ^ 1 - ^O, 


(4.21) 


(4.22) 


V {ran) ) 


fe=2 


where we have denoted by 2 (t) := 

/ {z{t),4i) \ 
{z{t),42) 

, t > 0 ;byA := 

/ Ai 0 . 

0 A 2 

• \ 
. 0 


V {z{t),4N) ) 


0 0 . 

• Aat / 


by(9:= 


/ {z{t),4i) \ 

0 


V 


0 


. Then, likewise in (14.1411 . we get 


/ 


-WA^Zml < -27iP^Z(t)||^ + S\\A\\\\Z{t)\\%, t > 0, 


(4.23) 


where ||A|| stands for the classical euclidean norm of the matrix A. Denote by 
Ai(^) > 0 the first eigenvalue of A, then by integration with respect to time in 
(I4.23L it yields 

\,{A)\\Zit)\\%<e-^^^^\\AiZo\\%+ re-2^^(‘-^)<5||Al||||Z(s)||^ds, (4.24) 

^0 

where making use of the Gronwall’s lemma 


\\mrN< 


1 


Ai(kl) 


\\A^Zo\\%exp 


(m\\ 

lUi(AI) 


t > 0 . 


(4.25) 


Let us denote by , i, j = 1,..., A^, the entries of the matrix i?i + 132 + ■•■ + Ba- 
By the definition of Bi (see (12.171) 1 and of the constants 7 ^, i = 1, ...,iV and S (see 
after (I2.13|) l. we have lim.y^. 
the adjoint of the matrix Bi 


^^^aoJi\bij\ G M+,Vi,j = 1,...,7V. Set 6 *- the entries of 


. + Bn. By definition of the adjoint and the above 
observation we deduce that lim-^j^oo | G M+, Vi,j = l,...,iV. Besides 


;and 






























STABILIZATION OF PARABOLIC SEMILINEAR EQUATIONS 


15 


this, the above observation also implies that lim^^ ->oo7i*' \det{Bi +... +Bpi)\ — 
+C)0. In other words, we have 


l^*, l < Ci, 


U 2(Af-l)’ 
7l 


b j = 1, and 


dct{B\ -j- 


B 


N) 


, 2(IV+1) 

< C7i 


for some positive constants Cij,c, = l,...,iV, independent of 71, for 71 large 
enough. This yields that, denoting by aij, i,j = 1 ,..., N, the entries of the matrix 
A = (Bi + ... + B]v)~^, there exists some constant C > 0 , independent of 71, such 
that \aij \ < i,j = 1 , ...,N. Consequently 


Pll < 


for 7i large enough. In conclusion, for 71 large enough, there exists some /i > 0 , 
such that 


Ai(yl) 


271 


ll^ll 

7 fAi(A) 


271 < CN ^ 


1 

A^ 


271 < -M, 


since xifA) —t 0 for 71 —>• 00. This, together with ( 14 . 251 ) . yields 


\\mrN< 


Ai(Al) 


,-Aii 






t > 0, 


that is the exponential decay of the first N modes of z. The rest of the proof follows 
immediately as the proof of Theorem 12.11 and it is omitted. □ 

Finally, for the proof of Theorem 12.21 similarly as for the one of Theorem 12.41 
one may proceed as in [U Theorem 4.2], or as in [TTl Theorem 1.1]. Therefore, they 
are omitted. 


5. Appendix 

Lemma 5.1. Under assumption (H), for any p < 71 < 72 < ... < 7jV) we have 


1_ 1 _1 


71 jAi 

72 —Ai 

71 Y ^2 

72 —A2 

7 iYAiv 
72—A^^ 

7^0. 

( 5 . 1 ) 

1 

7 iV—Ai 

1 

7 A'‘ — A2 

1 

7iV—Aiv 




Proof. Let us prove this by mathematical induction over N. Step 1 , for A = 2 , we 
have 


1 1 

71Y'^2 
72 —Ai 72 —A2 


_ (Ai - A 2 )( 7 i - 72) 

(71 - ^i )(72 - A 2)(72 - Ai)(7i 


A 2 ) 


^ 0 > 


since Ai < A 2 and 71 < 72 . 
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Step 2, we assume that for iV — 1 the claim is true and prove it for N. To this 
end we have 


1 

1 

1 

71^^! 

7i-_A2 

7i — Aiv 

72— Ai 

72 —A 2 

72 —Aiv 

1 

1 

1 

7iV-Ai 

7N‘ —-^2 

7iv — Aiv 


(Subtracting from the first column the one, from the {N — 1)*^ column the one ) 


A^^ —Ai 

Aat — A 2 

1 

1 

1 

( 71 —Ai)( 7 i—Aat) 
Aat —Ai 

(71 —A 2 )( 7 i —Aiv) 

Aat — A 2 

( 71 —Aai-i)(7i — Aa^ ) 

Aai—Aai -1 

71—Aa- 
1 

( 72 —Ai )(72 —Aat) 

(72 —A 2 )( 72 —Aiv) 

(72 —Aai-i)(72—Am-) 

72—Aa 

Aat —Ai 

Aat-A 2 

Aai—Aat-i 

1 

(7iV—A i)(7^^—A^^) 

(7iv —A 2 )( 7 n- —Aiv) 

(TN — A AT- 1 ) (7N-— An ) 

7A —Aa 


(-l)A-l 

7Ar — 7fc “ 


71-Ai 

71 -A 2 

1 

1 

1 

1 

1 

72 —Ai 

72 —A 2 

72 —Aa -1 

1 

1 

1 


JN—^l IN—^2 


lN — ^N-1 


1 

1 

1 


(Subtracting the line from the first one, the line from the [N — \y^ one) 

11 1 


(-1)^-^ T^ hAA-Afc)(7A-7fe) 
In — (7fc ~ ^n){'1n — Afc) 


1 

72—Ai 


71— A 2 
1 

72 — A 2 
1 


7iV-i —Ai 7iV-i —A 2 


71 —Aiv-i 

1 

72 —A^^_i 

1 

7iV-i—Aiv-: 


^0 


by the inductive hypothesis and the fact that {Ai} < A 2 < ... < Xn and 71 < 72 < 

... < 7Ar. □ 


Lemma 5.2. The sum Bi + B 2 + ... + Bn is an invertible matrix, where Bk, k = 
1, ...,N are introduced in relation \2.6\) . 


Proof. Arguing by contradiction, let us assume that there is z 
nonzero, such that (i?i + ... + Bn)z = 0. It follows that 

N 

{Bkz,z)j^ = 0 , 

fc=l 


Zi \ 

G K^, 


\ ZN J 


or, equivalently. 



Zi 


1 d 

7fe - Aj diy 



2 

dx = 0. 


Because the Lebesgue measure of Ti is nonzero, we deduce from the above that 


N 




1 d 

7 fe - dv 


(j),{x) = 0, 


a.e. on Ti, 


for all k = This gives N x N linear homogeneous systems, with the un¬ 

knowns Zi, i = 1 ,..., N, for almost all x G Ti. By the unique continuation property 
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of the eigenfunctions we know that, for all z = 1, 2, N, 7 ^ 0 on Fi. 

The determinant of the matrix of the corresponding system is 


1 

d 

7i-Ai 

dv' 

1 

d , 

72— 

dv' 

1 

d 


dv 


(j)l[x) 


N 


= na;^*( 


2=1 


1 1 


71 jAi 

7 l -^-^2 

to 

1 

72 — ^2 

1 

1 

7 N ‘ —-^1 

7 iV — A2 


1 


71 —Aiv 

12—^N 


IN — ^N 


/ 0 , a.e. X € Fi, 


by Lemma l5.1l Hence, necessarily z = 0. This is in contradiction with our assump¬ 
tion. We conclude that the sum Hi -h ... + is indeed an invertible matrix. □ 
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